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ABSTRACT 

We present a statistical study of the detectability of lightcurves of Kuiper Belt objects (KBOs). Some 
Kuiper Belt objects display lightcurves that appear "flat", i.e., there are no signiflcant brightness 
variations within the photometric uncertainties. Under the assumption that KBO lightcurves are 
mainly due to shape, the lack of brightness variations may be due to (1) the objects have very nearly 
spherical shapes, or (2) their rotation axes coincide with the line of sight. We investigate the relative 
importance of these two effects and relate it to the observed fraction of "flat" lightcurves. This study 
suggests that the fraction of KBOs with detectable brightness variations may provide clues about the 
shape distribution of these objects. Although the current database of rotational properties of KBOs is 
still insufScient to draw any statistically meaningful conclusions, we expect that, with a larger dataset, 
this method will provide a useful test for candidate KBO shape distributions. 
Subject headings: Kuiper Belt — minor planets, asteroids — solar system: general 



1. INTRODUCTION 

The Kuiper Belt holds a large population of small ob- 
jects which are thought to be r emnants of the proto- 
solar nebula ijJewitt fc Luu 1993|) . The Belt is also the 
most likely origin of other outer solar system objects such 
as Pluto-Charon, Triton, and the short-period comets; 
its study should therefore provide clues to the under- 
standing of the processes that shaped our solar system. 
More than 650 Kuiper Belt objects (KBOs) are known 
to date and a total of about 10^ objects larger than 
50 km are thought to orbit the Sun beyond Neptune 
llJewittfc Luu 20001). 

One of the most fundamental ways to study phys- 
ical properties of KBOs is through their lightcurves. 
Lightcurves show periodic brightness variations due to 
rotation, since, as the KBO rotates in space, its cross- 
section as projected in the plane of the sky will vary due 
to its non-spherical shape, resulting in periodic bright- 
ness variations (see Fig. A well-sampled lightcurve 
will thus yield the rotation period of the KBO, and 
the lightcurve amplitude has information on the KBO's 
shape. This technique is commonly used in planetary as- 
tronomy, and has been developed extensively for the pur- 
pose of determining the shapes, internal density struc- 
tures, rotational states, and surface properties of atmo- 
sphereless bodies. These properties in turn provide clues 
to their formation and coUisional environment. 

Although lightcurves studies have been carried out rou- 
tinely for asteroids and planetary satellites, the num- 
ber of KBO lightcurves is still meager, with few of suf- 
ficient quality for analysis (see Table ^1. This is due 
to the fact that most KBOs are faint objects, with ap- 
parent red magnitude of mR ~23 (Truiillo et al. 2001), 
rendering it very difficult to detect small amplitude 
changes in their brightness. One of the few high qual- 
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Fig. 1. — The lightcurve of an ellipsoidal KBO observed at aspect 
angle 9 = -k 12. Cross-sections and lightcurve are represented for 
one full rotation of the KBO. The amplitude. Am, of the lightcurve 
is determined for this particular case. See text for the general ex- 
pression. 



ity lightcurves is that of (20000) Varuna, which shows 
an amplitude of Am = 0.42 ±0.02 mag and a period 
of Prot = 6.3442 ± 0.0002 hrs l|Jewitt fc ShcDDard 2002). 
Only recently have surveys started to yield significant 
numbers of KBO s bright enough for detailed studies 
ijJewitt et al. 199 8). 

Another difficulty associated with the measurement of 
the amplitude of a lightcurve is the one of determining 
the period of the variation. If no periodicity is appar- 
ent in the data, any small variations in the brightness of 
an object must be due to noise. Furthermore, a precise 
measurement of the amplitude of the lightcurve requires 
a complete coverage of the rotational phase. Therefore, 
any conclusion based on amplitudes of lightcurves must 
assume that their periods have been determined and con- 
firmed by well sampled phase plots of the data. 

However, not all of the observed KBOs show detectable 
brightness variations (the so-called "flat" lightcurves). 
The simplest explanations for this could be due to (1) 
the object is axisymmetric (the two axes perpendicular 
to the spin vector are equal), or (2) its rotation axis is 
nearly coincident with the line of sight (see Fig. . In 
other words, the undetectable variations are either a con- 
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Table 1. KBOs with measured lightcurves. 



Name 


Class'' 


H° 




pa 


Source^ 






[mag] 


[mag] 


[hrs] 




1 qqq 


Q 


6.9 


<0.04 




RT99 


1994 TB 


p 


7.1 


0.3 


6.5 


RT99 


1996 TLee 


S 


5.4 


<0.06 




RT99 


1996 TPee 


P 


6.8 


<0.12 




RT99 


1994 VKg 


C 


7.0 


0.42 


9.0 


RT99 


1996 T066 


c 


4.5 


0.1 


6.25 


HaOO 


Varuna 


c 


3.7 


0.42 


6.34 


JS02 


1995 QYg 


p 


7.5 


0.6 


7.0 


RT99 


1996 RQ20 


c 


7.0 






RT99 


1996 TS66 


c 


6.4 


<0.16 




RT99 


1996 TQ66 


c 


7.0 


<0.22 




RT99 


1997 CS29 


c 


5.2 


<0.2 




RT99 


1999 TDio 


s 


8.8 


0.68 


5.8 


CoOO 



''dynamical class (C - classical KBO, P - plutino, S - scattered 
KBO) 
'^absolute magnitude 
"^liglitcurve amplitude 
'^rotatio nal period 

^TS02 I.Tewitt fc Sheppa-rd 20021) HaOO IHainaut 200(ll. C oOO 
IConsolmagno et al. 20001) . RT99 IRomanishin Tegler 19931 



Table 2. Used symbols and notation. 



Symbol Description 

a>b>c axes of ellipsoidal KBO 

a>b>c normalized axes of KBO (6 = 1) 

aspect angle 

Arrtniin minimum detectable lightcurve amplitude 

Bmin aspect angle at which Am = Ammin 



sequence of the KBO's shape, or of the observational ge- 
ometry. By studying the relative probabilities of these 
two causes, and relating them to the observed fraction of 
"flat" lightcurves, we might expect to improve our knowl- 
edge of the intrinsic shape distribution of KBOs. In this 
paper we address the following question: Can we learn 
something about the shape distribution of KBOs from 
the fraction of "flat" lightcurves? 

2. DEFINITIONS AND ASSUMPTIONS 

The observed brightness variations in KBO lightcurves 
can be due to: 

• eclipsing binary KBOs 

• surface albedo variations 

• irregular shape 

In general the brightness variations will arise from some 
combination of these three factors, but the preponder- 
ance of each effect among KBOs is still not known. In 
the following calculations we exclude the first two factors 
and assume that shape is the sole origin of KBO bright- 
ness variations. We further assume that KBO shapes 
can be approximated by triaxial ellipsoids, and thus ex- 
pect a typical KBO lightcurve to show a set of 2 maxima 
and 2 minima for each full rotation (see Fig.nj. Table El 
summarizes the used symbols and notations. The listed 
quantities are defined in the text. 
The detailed assumptions of our model are as follows: 



a) 



line of sight 




Fig. 2. — a) A spherical coordinate system is used to represent 
the observing geometry. The line of sight (oriented from the object 
to the observer) is the polar axis and the azimuthal axis is arbitrary 
in the plane orthogonal to the polar axis. 9 and </) are the spherical 
angular coordinates of the spin axis s. In this coordinate system 
the aspect angle is given by 9. The "non-detectability" cone, with 
semi- vertical angle Smim is represented in grey. If the spin axis lies 
within this cone the brightness variations due to changing cross- 
section will be smaller than photometric errors, rendering it im- 
possible to detect brightness variations, b) The picture represents 
an ellipsoidal KBO with axes a > b > c. 



1. The KBO shape is a triaxial ellipsoid. This is the 
shape assu med by a rotating bod y in hydrostatic 
equilibrium ijChandrase khar 1969). There are rea- 
sons to believe tha t KBOs might have a "rub- 
ble pile" structure l|Farinella et al. 19811) . justify- 
ing the approximation even further. 

2. The albedo is constant over surface. Although 
albedo variegation can in p rinciple explain any 
given lightcurve tR ussell 1 9Q6[), the large scale 
brightness variations are generally attributed to t he 
object's irregular shape l|Burns fc Tedesco 19791) . 

3. All axis orientations are equally probable. Given 
that we have no knowledge of preferred spin vec- 
tor orientation, this is the most reasonable a priori 
assumption. 

4. The KBO is in a state of simple rotation around 
the shortest axis ( the axis of maximum moment of 
inertia). This is likely since the damping timescale 
of a complex rotation (e .g. , precession) , ^ 10"^ yr, 
ijBurns fc Safronov 1973(1 . l|Harris 1994 is smaller 
than the estimated time between collisions (lO'^- 
10"'^^ yr) that would re-excite such a rotational state 
fStern 1995), (Davis & Farinella 1997). 

5. The KBO is observed at zero phase angle {a — 0) It 
has been shown from asteroid data that lightcurve 
amplitudes seem to increase linearly with phase an- 
gle, 

A{e,Q) = A{6,a)/{l+ma) , 

where 9 is the aspect angle, a is the phase angle and 
TO is a coefficient which depends on surface compo- 
sition. The aspect angle is defined as the angle 
between the line of sight and the spin axis of the 
KBO (see Fig. |2K) , and the phase angle is the Sun- 
object-Earth angle. The mean values of m found 
for different asteroid classes are to(S) = 0.030, 
to(C) = 0.015, m{M) = 0.013, where S, C , and M 
are asteroid classes ((Michalowskv 199^ . Since 
KBO are distant objects the phase angle will al- 
ways be small. Even allowing m to be one order 
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of magnitude higher than that of asteroids the in- 
crease in the hghtcurve amphtude wiU not exceed 
1%. 

6. The brightness of the KBO is proportional to its 
cross-section area (geometric scattering law). This 
is a good approximation for KBOs because (1) most 
KBOs are too small to hold an atmosphere, and (2) 
the fact that they are observed at very small phase 
angles reduces the influence of scattering on the 
lightcurve amplitude (Magnusson 1989). 

The KBOs will be represented by triaxial ellipsoids 
of axes a > b > c rotating around the short axis c (see 
Fig-Eb)- In order to avoid any scaling factors we normal- 
ize_all axes by b, thus obtaining a new set of parameters 
a, b and c given by 



a — a/b , b — 1 



c/b 



(1) 



As defined, a and c can assume values 1 < a < cx) and 
< c < 1. Note that the parameters a and c are dimcn- 
sionless. 

The orientation of the spin axis of the KBO relative to 
the line of sight will be defined in spherical coordinates 
{0,4>), with the line of sight (oriented from the object 
to the observer) being the z-axis, or polar axis, and the 
angle 6 being the polar angle (see Fig. [5^). The solution 
is independent of the azimuthal angle (/), which would 
be measured in the plane perpendicular to the line of 
sight, between an arbitrary direction and the projection 
of the spin axis on the same plane. The observation 
geometry is parameterized by the aspect angle, which in 
this coordinate system corresponds to 0. 

As the object rotates, its cross-section area S will vary 
periodically between S'max and 5111111 (see Fig. Pl. These 
areas are simply a function of a, &, c and the aspect angle 
0. Given the assumption of geometric scattering, the 
ratio between maximum and minimum flux of reflected 
sunlight will be equal to the ratio between 5niax and S'min- 
The lightcurve amplitude can then be calculated from the 
quantities a, c and 9 and is given by 



Am = 2.5 log 



cos^ 9 + a'^c'^ sin^ ( 
52 cos^ 9 + (? sin^ 9 



1/2 



(2) 



3. "FLAT" LIGHTCURVES 



It is clear from Eq. |5J that under certain condi- 
tions. Am will be zero, i.e., the KBO will exhibit a flat 
lightcurve. These special conditions involve the shape 
of the object and the observation geometry, and are de- 
scribed quantitatively below. Taking into account pho- 
tometric error bars will bring this "flatness" threshold to 
a finite value, Am,nin, a minimum detectable amplitude 
below which brightness variation cannot be ascertained. 

The two factors that influence the amplitude of a KBO 
lightcurve are: 

1. Sphericity For a given ellipsoidal KBO of axes 
ratios a and c the lightcurve amplitude will be largest 
when 9 = tt/2 and smallest when 9 = or tt. At 9 = 7r/2, 
Eq. becomes 



Am = 2.5 log a . 



(3) 



Even at 9 — 7r/2, having a minimum detectable ampli- 
tude, Amniin, puts constraints on a since if a is too small, 




Fig. 3. — Illustration of a rotating ellipsoid at different aspect 
angles. A quarter of a full rotation is represented. Rotational phase 
of ellipsoid increasing from top to bottom and 6 decreasing from 
left to right. T is the period of rotation. Axes ratios are a = 1.2 
and c = 0.9. 



the lightcurve amplitude will not be detected. This con- 
straint is thus 



d < 10 



0.4Am„ 



"flat" lightcurve . 



(4) 



2. Observation geometry If the rotation axis is 
nearly aligned with the line of sight, i.e., if the aspect 
angle is sufficiently small, the object's projected cross- 
section will hardly change with rotation, yielding no de- 
tectable brightness variations (see Fig. The finite 
accuracy of the photometry defines a minimum aspect 
angle, ^min: within which the lightcurve will appear flat 
within the uncertainties. This angle rotated around the 
line of sight generates the "non-detectability cone" (see 
Fig.l^t), with the solid angle 



^^(^min) 



sin 9 d9 ( 



(5) 



Any aspect angle 9 which satisfies 9 < ^min falls 
within the "non-detectability cone" and results in a non- 
detectable lightcurve amplitude. Therefore, the proba- 
bility that the lightcurve will be flat due to observing 
geometry is 

2 X n{9^in) 



Pa c (non-detection) : 

= 1 — cos 6 
Pa,E(detection) =cos 

"mill- 



(6a) 
(6b) 

The factor of 2 accounts for the fact that the axis might 
be pointing towards or away from the observer and still 
give rise to the same observations, and the 4tt in the 
denominator represents all possible axis orientations. 

From Eq. ^ we can write cos ^min as a function of a 
and c, 



COS! 



*(a, c) 



c2(a2 _ K) 



(7) 



where K = iQO S^m^^in^ xhe function *(a, c), repre- 
sented in Fig. 21 is the probability of detecting bright- 
ness variation from a given ellipsoid of axes ratios (a, c). 
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+ (a,c) 




Fig. 4.— The function 1'(a, c) (Eq. Q). This plot assumes 
photometric errors Am^i^ = 0.15 mag. The detection probabihty 
is zero when a < loO-^^^min ~ 1.15. 



It is a geometry weighting function. For a in [1, v^^] we 
have ^(a, c) = by definition, since in this case the KBO 
satisfies Eq. Q and its hghtcurve amphtude will not be 
detected irrespective of the aspect angle. It is clear from 
Fig.Elthat it is more likely to detect brightness variation 
from an elongated body. 

4. DETECTABILITY OF LIGHTCURVES 

In order to generate a "non-flat" lightcurve, the KBO 
has to satisfy both the shape and observation geometry 
conditions. Mathematically this means that the proba- 
bility of detecting brightness variation from a KBO is a 
function of the probabilities of the KBO satisfying both 
the sphericity and observing geometry conditions. 

We will assume that it is possible to represent the 
shape distribution of KBOs by two independent prob- 
ability density functions, /(a) and g{c), defined as 

1*0,2 roo 

p(ai<a<a2)= / /(a)da, / /(a)da=l, (8a) 

Jai Jl 

p{ci<c<C2)^ r g{c)dcj <?(c)dc=l, (8b) 

Jcx Jo 

where the integrals on the left represent the fraction of 
KBOs in the given ranges of axes ratios. This allows us 
to write the following expression for p(Am > Amiuin), 
where both the shape and observation geometry con- 
straints are taken into account. 



p{Am > Aminin) 

poo 







^'(a, c)f{a)g{c) da dc. 



(9) 



The right hand side of this equation represents the proba- 
bility of observing a given KBO with axes ratios between 
(a, c) and (a -I- da, c -I- dc), at a large enough aspect an- 
gle, integrated for all possible axes ratios. This is also 
the probability of detecting brightness variation for an 
observed KBO. 
The lower limit of integration for a in Eq. can be 

replaced by Vk, with K defined as in Eq. {Tj), since 
\l/(a, c) is zero for a in [1, ^/K]. In fact, this is how the 
sphericity constraint is taken into account. 



Provided that we know the value of p(Am > Ammin) 
Eq. © can test candidate distributions /(a) and g{c) 
for the shape distribution of KBOs. The best estimate 
for p(Am > Amniin) is given by the ratio of "non- 
flat" lightcurves {N-q) to the total number of measured 
lightcurves (TV), i.e., 



p(Ato > Ammiii) 



N ' 



(10) 



Because N is not the total number of KBOs there will 
be an error associated with this estimate. Since we do 
not know the distributions /(a) and g{c) we will assume 
that the outcome of an observation can be described by 
a binomial distribution of probability p(Am > Arnmin). 
This is a good approximation given that N is very small 
compared with the total number of KBOs. Strictly 
speaking, the hypergeometric distribution should be used 
since we will not unintentionally observe the same ob- 
ject more than once (sampling without replacement). 
However, since the total number of KBOs (which is not 
known with certainty) is much larger than any sample of 
lightcurves, any effects of repeated sampling will be neg- 
ligible, thereby justifying the binomial approximation. 
This simplification allows us to calculate the upper (p+) 
and lower (p_) limits for p(Am > Am,„iii) at any given 
confidence level, C. These values, known as the Clopper- 
Pearson confidence limits, can be f ound solving t he fol- 
lowing equations by trial and error l)Barlow 19891) . 



N 



P{r;p+{Am>Am,^in),N) 



r=ND + l 



J2 P{r;p-iAm>Am^in),N) 



C + 1 



C+1 



(Ua) 



(lib) 



r=0 



(see Table |5| for notation) where C is the desired con- 
fidence level and P{r;p,N) is the binomial probability 
of detecting r lightcurves out of N observations, each 
lightcurve having a detection probability p. Using the 
values in Tableland Ammin — 0.15 mag we have Nd = 5 
and iV = 13 which yields p{Am > Am^in) = 0.38t°-J^ 
at a C = 0.68 (Icr) confidence level. At C = 0.997 
(3cr) we have p{Am > Ammin) — 0.38lo.3i- The value 
of p{Am > Amniin) could be smaller since some of the 
flat lightcurves might not have been published. 

Note that for moderately elongated ellipsoids (small 
a) the function ^'(a, c) is almost insensitive to the pa- 
rameter c (see Fig. , in which case the axisymmetric 
approximation with respect to a can be made yielding 
c « 1. Eq. then has only one unknown parameter, 
/(«)■ 



p{Am > Ammin) ' 



0.38 



«'(a,l)/(a)da 



-1-0.41 
0.31- 



(12) 



If we assume the function /(a) to be gaussian, we can use 
Eq. (|12|l to determine its mean /z and standard deviation 
(7, after proper normalization to satisfy Eq. JHaJ. The 
result is represented in Fig. |S1 where we show all pos- 
sible pairs of (/i,cr) that would satisfy a given p(Am > 
Ammin)- For example, the line labeled "0.38" identi- 
fies all possible pairs of (/i,cr) that give rise to p{Am > 
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Theoretical _p(Am>Am^[^) 
1 .40 F ' ' ' ' ' ' ' ' "x^ 




Fig. 5. — Contour plot of the theoretical probabilities of detect- 
ing brightness variation in KBOs (assuming Am^^^ = 0.15 mag), 
drawn from gaussian shape distributions parameterized by fi and a 
(respectively the mean and spread of the distributions). The solid 
lines represent the observed ratio of "non-fiat" lightcurves (at 0.38) 
and 0.68 confidence limits (at 0.23 and 0.56 respectively). 



Non-uniform albedo would cause nearly spherical 
KBOs to generate detectable brightness variations, de- 
pending on the coordinates of the albedo patches on the 
KBO's surface. This means that our method would over- 
estimate the number of elongated objects by attributing 
all brightness fluctuations to asphericity. 

Binary KBOs would influence the results in different 
ways depending on the orientation of the binary system's 
orbital plane, on the size ratio of the components, and 
on the individual shapes and spin axis orientations of 
the primary and secondary. For example, an elongated 
KBO observed equator-on would have its lightcurve flat- 
tened by a nearly spherical moon orbiting in the plane 
of the sky, whereas two spherical KBOs orbiting each 
other would generate a lightcurve if the binary would be 
observed edge-on. 

These effects are not straightforward to quantify an- 
alytically and might require a different approach. We 
intend to incorporate them in a future study. Also, with 
a larger sample of lightcurves it would be useful to apply 
this model to subgroups of KBOs based on dynamics, 
size, etc. 



Ammin) — 0.38, the line labeled "0.56" all possible pairs 
of (/x,cr) that give rise to p{Am > Ammin) = 0.56, etc. 

Clearly, with the present number of lightcurves the 
uncertainties are too large to draw any relevant conclu- 
sions on the shape distribution of KBOs. With a larger 
dataset, this formulation will allow us to compare the 
distribution of KBO shapes with that of the main belt 
asteroids. The latter has been shown to resemble, to 
some extent, that of fragments of high-velocity impacts 
ijCatullo et al. 198^ . It deviates at large asteroid sizes 
that have presumably relaxed to equilibrium figures. A 
comparison of /(a) with asteroidal shapes should tell us, 
at the very least, whether KBO shapes are collisionally 
derived, as opposed to being accretional products. 

The usefulness of this method is that, with more data, 
it would allow us to derive such quantitative parameters 
as the mean and standard deviation of the KBO shape 
distribution, if we assume a priori some intrinsic form 
for this distribution. The method's strength is that it 
relies solely on the detectability of lightcurve amplitudes, 
which is more robust than other lightcurve parameters. 

This paper focuses on the influence of the observation 
geometry and KBO shapes in the results of lightcurve 
measurements. In which direction would our conclusions 
change with the inclusion of albedo variegation and/or 
binary KBOs? 
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5. CONCLUSIONS 

We derived an expression for the probability of detect- 
ing brightness variations from an ellipsoidal KBO, as a 
function of its shape and minimum detectable amplitude. 
This expression takes into account the probability that 
a "flat" lightcurve is caused by observing geometry. 

Our model can yield such quantitative parameters as 
the mean and standard deviation of the KBO shape dis- 
tribution, if we assume a priori an intrinsic form for this 
distribution. It concerns solely the statistical probabil- 
ity of detecting brightness variation from objects drawn 
from these distributions, given a minimum detectable 
lightcurve amplitude. The method relies on the assump- 
tion that albedo variegation and eclipsing binaries play a 
secondary role in the detection of KBO lightcurves. The 
effect of disregarding albedo variegation in our model is 
that we might overestimate the fraction of elongated ob- 
jects. Binaries in turn could influence the result in both 
directions depending on the geometry of the problem, 
and on the physical properties of the constituents. We 
intend to incorporate these effects in a future, more de- 
tailed study. 
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